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Abstract 

The energy dependence of the nonradiative electron capture cross- 
section is discussed in the relativistic domain. A simple analytic ex- 
pression is obtained for inner-shell transitions using second-order per- 
turbation theory. We have confirmed that the leading-order term is 
found to have the following energy dependence: a ~ E- X \\^E. This 
is attributed to a combination of kinematic features of the process 
and retardation effects. Electron capture without change of spin is 
the dominant transition. 

PACS 31.15.Md, 31.30.Jv, 34.70+e, 52.20.Hv 

The asymptotic form of the electron capture process for relativistic ener- 
gies has been a subject of debate for some time. An overview of this ques- 
tion has been given by Mukherjee et al (1992) and Bransden and McDowell 
(1993), while a detailed discussion is given by Eichler and Meyer hof (1995). 
In the relativistic domain, the physical constraints of energy and momentum 
matching conspire to make capture a weakly-coupled process in comparison 
to ionization or excitation. Indeed it is interesting to remark that for rela- 
tivistic energies, and for small nuclear charges, vacuum coupling in the form 
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of radiative electron capture is the dominant mechanism for charge exchange. 
As a consequence, the numerical results are primarily of academic interest 
in the extreme relativistic domain. Nonetheless, results for the nonradiative 
process ultimately provide important benchmarks and calibration standards 
for both experiment and theory (Eichler and Meyerhof 1995). Moreover the 
phenomenon has attracted attention in recent years not least because it is 
physically important for large nuclear charges, and this makes it accesssible 
to experimental studies up to and beyond the energy range GeV/u (Anholt 
1985). In the following discussion we use the Lorentz factor 7 as the energy 
parameter, and this is related to the laboratory frame kinetic energy through 
the usual relation E — (7 — l)Mpc 2 , with Mp the projectile mass. 

Experimental data for relativistic collisions have been analysed and com- 
pared with theoretical models (Eichler 1990, Moiseiwitsch 1989, Deco and 
Rivarola 1987, Glass et al 1994), and satisfactory agreement has been ob- 
tained in many cases. However an intriguing question still remains open 
regarding the high-energy limit of the cross section, both in terms of the en- 
ergy dependence and the size of the cross section. The uncertainty over this 
result reflects a well-known anomaly for high-energy nonrelativistic collisions, 
namely that electron capture is dominated by second-order terms. Nonethe- 
less the dominance of the second-order term, via the Thomas mechanism, 
does not come into effect until the energies enter the relativistic domain. 
Paradoxically, as the energy rises further this mechanism becomes less im- 
portant than the first-order process due to the contraction of suitable phase 
space. In fact it has been shown that the contribution to the cross section 
(a) given by the Thomas term decreases as a ~ 7 -3 (Shakeshaft and Spruch 
1979). This is much faster than the first-order term which behaves as a ~ 7 _1 
(Shakeshaft 1979) and ultimately it becomes negligible in the ultrarelativistic 
region. 

In applying the Bates refinement of the ROBK1 approximation (Moisei- 
witsch and Stockman 1980), the results were found to give rather good cross 
section estimates for the intermediate energy range. However this model also 
led to the prediction that a ~ 7~ 1 (hi7) 2 at extremely high energies. The 
result that the nonorthogonal correction term eventually dominates the pri- 
mary term (ROBK1) was suprising, and subsequently misgivings have been 
expressed regarding the validity of this approach (Eichler 1990). This result 
also contradicted other modified first-order models, such as relativistic BIB 
(Moiseiwitsch 1986) which predicted the result a ~ 7 _1 (Toshima and Eich- 
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ler 1990). Conversely it could be argued that the important influence of this 
correction term is not more suprising than the dominance of second-order 
terms, and reflects the weakness of the first-order process. 

The picture in second-order theories is not so clear. Most of the work 
to date has centred on the well-known ROBK2 model (Moiseiwitsch 1985, 
Decker 1990) and relativistic impulse approximation RIA (Jakubassa- Amundsen 
and Amundsen 1980). Humphries and Moiseiwitsch (1985) followed the stan- 
dard nonrelativistic approach to evaluating ROBK2 by using a linearised- 
propagator peaking approximation as well as making expansions in aZ in 
the wave functions. This approximation reduces the T-matrix to a simple 
analytic form, and gives the result a ~ as in ROBK1 though with a 
modified coefficient. Decker and Eichler (1993) managed to calculate the T- 
matrix exactly using numerical techniques, and in doing so queried the earlier 
ROBK2 work on two counts. Firstly, when they compared their results with 
experimental data (which involves largish charges) they found that ROBK2 
was in worse accord with experiment than ROBK1. This discrepancy has 
not as yet been fully accounted for in that the relative effects of the aZ ex- 
pansions and the peaking approximation have not been isolated. Decker and 
Eichler (1993) also claimed signs of an energy dependence which contained 
logarithmic terms. Again possibly indicative of the size of the correspond- 
ing coefficient, the evidence was not conclusive even though they performed 
calculations for 7 > 1000. 

Jakubassa- Amundsen and Amundsen (1980) in extending their work on 
the impulse approximation to the relativistic regime had argued that the 
ionization-like matrix element embedded in the RIA T-matrix would lead to 
logarithmic terms. The origin of these terms is the retardated potential. In- 
deed logarithmic correction terms are ubiquitous in high-energy atomic pro- 
cesses such as the stopping of relativistic charged particles. In both ROBK2 
and RIA the possibility of retardation is obvious from the form of the T- 
matrix but as yet there has been no detailed mathematical study of both the 
origin of these terms as well their size in comparison to the ROBK1 ampli- 
tude. We have attempted to this effect. We will show that the multiplicative 
factor upon the logarithmic terms is of order a 2 (Zp + Zt) 2 compared with 
ROBK1. 
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1 First-order Approximations 



The cross section for electron capture (a), for a given impact energy in the 
laboratory frame, can be written in terms of the scattering amplitude T(rj) 
as follows (Glass et al 1996): 

(7(7) = (^IvY 2 j 'dr, \T(rj)\ 2 . (1) 

In this formula, v is the relative speed of target (T) and projectile (P), and 
7 = (1 — t> 2 /c 2 ) -1 / 2 . The momentum transfer variable is denoted by t?, and 
all quantities are expressed in atomic units (|e| = m e — h — 1, and a = 1/c). 

Specifically, we consider the capture of an electron from the target ls- 
state $t to the projectile ls-state ($p). In this paper we use the following 
notation to describe these states: 

$ r = (t> T ( r )e- iETt , & P = <j>' P {r')e- iEp1f . (2) 

In the equation above, the functions are defined in the rest frames of the 
target and projectile nuclei respectively. The electron coordinate is denoted 
by (r) with respect to the origin located at the target nucleus, and it is 
labelled (r') when referred to the projectile nucleus. These coordinates are 
related through the usual Lorentz transformation. All primed variables in- 
dicate that they are associated with the rest frame of the projectile nucleus. 
The atomic energies, which include the electron rest mass, are denoted by 
E PT . Then the relativistic first-order OBK (ROBK1) approximation takes 
the form (Shakeshaft 1979): 

TrobkiM = (2tt) 3 / 2 J d 3 s 4>p(K' + s)V P (s)SMK). (3) 

The Lorentz transform operator for a boost from the target to the projectile 
frame is denoted by S. The tilde signifies the Fourier transform through the 
definition 

f(p) = (2n)-^Jd 3 re-V r f(r). (4) 
The momentum transfer vectors are simply given by: 

K = K ± + K\\ = 77 + v^Et - E P )/(iv), 

K' = K'j_ + K'u = T) — v(^E P — Et)/(^v). (5) 
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Evidently at very high velocities (7 — > 00) the only energy dependence 
in the expression for Trobki is due to the operator S, and thus T ~ 7 1 / 2 . 
Then it clear that in the same limit a ~ 7 _1 . This energy dependence is a 
feature of the kinematics and is due to a combination of time-dilation and 
frame transformation. The cross section for ground-state capture, without 
change of spin (orqbki) then takes the simple form, to leading order in aZ: 

-£>b K1 (7) « (6) 

in atomic units. 

The difficulty with calculating electron capture using this model is that 
the result, though Lorentz invariant, is not gauge invariant due to the nonortho£ 
onality of the intial and final states. Moiseiwitsch and Stockman (1980) 
investigated this problem using a modification of the ROBK1 originally 
due to Bates. This defines a correction due to orthogonalization given by: 
Tbrb = Trobki — T corr , where, 

Tcorriv) = (2vr) 3 J d 3 s $,{K' + s)S- 1 F(s)S 2 V^s)4> T (K + a). (7) 

In this expression 

F(s) = (2tt)- 3 / 2 J d 3 r e^ r 4(r)0 r (r). (8) 

We use the overline symbol to denote a stretched vector thus: s = s± + 
^SzV. In physical terms this correction term allows for capture through 
elastic scattering followed by the sudden adjustment of the wavefunction 
from target to projectile. Although this represents a two-step process, the 
matrix element is of first order in the interaction potential. In mathematical 
terms the new feature introduced in equation (^) is given by F(s) which is 
strongly peaked around s = and contains 7-dependent terms. The integral 
can be evaluated by the peaking approximation. One can show (Glass 1994) 
that the expression (|7|) is strongly determined by this peak, and to a high 
degree of accuracy one can simply replace the momentum distributions by 
their values at this peak: $p(K' + s) w (fip(K'), 4>t{K + a) w 4>t(K). 
Therefore the contribution of this term has a weight determined by the tails 
of the momentum distributions of the atoms, and as a result will be of higher 
order in aZ than the term Trobki- After making these assumptions, the 



5 



integral poses no further difficulty and we obtain the following result for 
capture without spin flip: 




Jco„(i7)~ ( v /c)(IC + Zl )''(K" + Zj,y \ 2c')- [ > 

The (In 7) term is a consequence of the the binary scattering peak V P (s) and 
the presence of this factor will ultimately dominate the energy dependence 
so that as 7 — > 00: 



*B»<7> « (^f) ^(Mf (10) 



in agreement with Moiseiwitsch and Stockman (1980) for capture from T(ls) 
to P(ls) without change of spin, and to leading order in aZp t T- As mentioned 
by these authors, such a term requires extremely high energies and/or large 
charges for its presence to be significant. 



2 Second-order corrections 

The second-order term of the ROBK perturbation series gives the result 
(Glass et al 1996): 

T 2 (rj) = J d 3 s J d 3 l 4>1{K' + l)V P (l)S- 1 
X — — 1,,. = — — ; — s 2 v t (s)4> t (k + S). (11) 



(3r + -iv ■ If - c 2 (K + I + s) 2 - c 4 + le 

defined in the limit e — > 0+. The Dirac matrices are denoted by ex. and /3. If a 
direct comparison is to be made with experimental data this T-matrix should 
be calculated without further approximations. However, in order to make 
analytic progress, and if we are primarily concerned with the asymptotic form 
of the cross section, we can take aZ to be a good perturbative parameter, 
and seek the leading-order term of this integral. Historically, simplifying 
assumptions have been used in order to reduce this integral to a simple 
analytic form (Bransden and McDowell 1993). In general these methods 
relied upon peaking approximations which are known to have limited scope. 
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As in section [I] the basic idea is that if an integrand possesses sharp and 
isolated peaks then the integral may be simplified. In expression (|TT|) there 
are four significant peaks in the integrand at the values: l+K' = 0, s+K = 0, 
and I = 0, s = 0. The first pair corresponds to bound-state momentum 
distributions while the second pair is associated with binary scattering. The 
fact that the peaks are sharp allows one to stretch the integration limits 
from the neighbourhood of the peak to infinity. And the isolation of the 
peaks permits one to separate each term. The leading term arises from the 
double bound-state peaks (I + K' — 0, s + K = 0). However there will be 
cross terms from the binary and bound-state peaks at (l = 0,s + K = t = 
) and at ( I + K' = T = 0, s = 0). Although these terms are are of higher 
order in powers of aZ, we find that that they have a slower decay as a 
function of increasing energy and therefore dominate in the ultrarelativistic 
limit. Finally there will be a term of still higher order in aZ that arises from 
the double binary peak (Z = 0, s = 0). On the understandng that the peaks 
are sharp and isolated, then the expresssion (|ll]) can be divided into four 
components as follows: 

T% = Thm + Tb pa + Tbpb + Tqbp (12) 

The first term (Thm) was evaluated by Humphries and Moiseiwitsch (1984). 
Retaining zeroth and first-order terms in I + K' and s + K in the propagator, 
they obtained the expression: 

^ 2^{Z P Z T fl 2 J p {-K)S-\-ccl ■ K + (3c 2 + 7 e / )^V T (K / ) 
hmW ~ K2RI2 ( 2l E T E P -c 4 -e}- K 2 c 2 + 2iZ P Kc 2 + 2iZ T K'c 2 ) 

(13) 

where cr(k) denotes a Darwin spinor 

CT TiP (k) = M + ^f^J VT,P. (14) 

In this notation the spin vector is defined as wp,T =(1 0) for spin up. 

For aZp t T C 1 we get 

T (n) - 2^{Z P Z T f/\^fl 2 F 
hmkJi) K 2 K ,2 { 2c 2 ( 1 -l)-K 2 + 2iZ P K + 2iZ T K') 1 ' 
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where F = 7 + | (^+f) f° r nonflip (TT) an d F = x{2^ + l)i]/2c for capture 
with spin flip (Ti)- 

This approximate form for the ROBK2 transition amplitude has to be 
qualified on several counts. Firstly, we would expect the peaking approxi- 
mations used in this derivation to break down for large charges and/or low 
velocities; neither criterion is strongly satisfied for much of the experimental 
data (Anholt 1985). However, it transpires that total cross sections within 
this approximation are in considerably better agreement with experiment 
than the ROBK1 results, which are commonly five times too large. On 
average the Humphries and Moiseiwitsch (HM) approximation to ROBK2 
is around a factor of two lower than ROBK1. This still leaves a consider- 
able gap between theory and experiment. Moiseiwitsch (1988) introduced 
the 'averaging approximation' in which the external K' 2 K 2 factors in the 
denominator of equation ( 15]) are replaced by K' 2 + Z\ and K 2 + Zj<, which 
resulted in a dramatic improvement between theory and experiment. 

In 1990 Decker calculated ROBK2 without approximation and using 
fully relativistic wavef unctions. It was discovered that in many cases rather 
than improving on the ROBK1 predictions, exact ROBK2 was an order of 
magnitude greater than ROBK1 (Eichler 1990), and consequently in poorer 
accord with experiment. This obviously raises serious doubts concerning the 
validity of the various peaking approximations discussed above, especially 
for large charges and energies in the lower range of the relativistic regime. 
There has been considerable debate in the literature as to which method is 
correct, with Decker and Eichler (1993) arguing in favour of exact numerical 
calculation of the perturbation series, while Moiseiwitsch (1995) has pointed 
out the possible pitfalls of simultaneously expanding along the Born series 
and in terms of aZ in the wavef unctions. 



2.1 Binary scattering peaks 

In section [I] it was shown that the binary scattering peak can produce loga- 
rithmic terms though with a small coefficent. The contribution of such terms 
within second-order perturbation theory can be evaluated without much dif- 
ficulty so long as the assumption that the peaks are sharp and isolated holds 
true. This supposes, for example, that the peaks at I = and I = —K' are 
distinct, which is certainly the case so long as aZp <C 1. In this calculation 
we continue to make this assumption, using aZ as a perturbation parameter, 
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in order to obtain the leading-order correction from the binary peaks. 
The first binary peak of (12) leads to the term: 

T BPA (r7) w [dH [ dH 4>*(K' + l)V P (l)S- 1 (16) 



x 



col ■ I + (3c 2 + E T + ■ I 



E%-c* + 2-/E T v ■ I - 2cH ■ I - c 2 l 2 + ie 



x S 2 V T (-K)4> T (t). 



and a similar term will aries from the second peak at: / + K' = T = 0, s = 0. 

In the usual manner, we can invert the Fourier transform for <f>r{t) by 
using the identity 



(A + ie) 



—i 



du exp(iu(X + ie)) 



(17) 



to convert the denominator into an exponential. Explicitly, we have for ls- 
states the expression: 



T BPA (*7) w (2^ 2 (Z" T /n) 1 / 2 ^(K')V T (-K) 
cot l + 13c 2 + E T + jv ■ I 



x 



E\ - c 4 + 2^E T v ■ I - c 2 l 2 + 2ic 2 Z T l + ie 



dH v' P (i)S- 1 

S 2 (T T (-l) 



To leading order in a Z we have 

T B pa(t7) = (2n) 3 / 2 (Zy7c) 1 / 2 ^(K')V T (-K)S- 1 ((3+l)S 2 w T L 

with 



J dH v' P {i)- 



l 



;is) 



(19) 



(20) 



(I 2 - 2v ■ I + Z\ - 2iZ T l - ie) 

Even after these simplifications, the calculation of this integral is not straight- 
forward. However, it is clear that the main feature of the integrand is due 
to the singularity in the propagator. It is this factor rather than the po- 
tential scattering in ( 20 ) that gives rise to the logarithmic behaviour. An 
analytic form can be found for this integral (Glass 1994) though its form is 
rather complicated. Instead it can be shown to a fair approximation that 
this complicated form reduces to, 



-i(2nf /2 (Z P /jv) ln(2ijZ T /v). 



(21) 
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Then within these approximations we have a simple formula for ground- 
state to ground-state capture without change of spin: 



7ffi A « -(2iZ P /v) H2irfv/Z T )T^ OBK1 . (22) 

Using the same method we find that the term due to the peaks around 
I + K' = T = 0, s = 0. is approximately given by: 

Tfia « -(JMZ^/t;) ln(2^/Z P )T^ BK1 . (23) 

While these terms are of higher order than ROBK1 in aZ, they are are 
of lower order than the Bates term given by equation (^) . Furthermore it is 
clear that in the limit of extreme relativistic energies that this term decreases 
much more slowly than the bound-state peaking term fllS]). If we combine 
these terms according to fll"2]) the result is that at relativistic energies the 
nonradiative electron capture cross section has the asymptotic form: 

TT M 2 7 irZS T Z P (Z T + Z P ) 2 (\n 1 y 
r ROBK 2 (7) ^ — (24) 

It is clear from this expression that the ROBK2 approximation contains 
logarithmic terms and that these are in general more important than the term 
from the Bates approximation (|H]) for small aZ. The simple reason is that 
only one momentum tail is sampled in producing the term (|T9|) . To complete 
our discussion, we note that additional logarithmic terms will be produced by 
the product of the two binary peaks: term Tdbp in equation (0). However, 
these will be of higher order in aZ than the term (O) because they include 
terms from both momentum tails as in the Bates approximation. Hence 
their contribution will not be a significant factor in this process for aZ <C 
1. Finally we conclude by mentioning that capture with spin-flip through 
this process contributes a small additional term. Further consideration of 
equation (|T9| ) for the case of capture with spin-flip shows that, in common 
with ROBK1, the cross section is only a quarter of the non-flip result (a^ = 
t(7^). Previously, calculations using ROBK2 (with bound-state peaking) 
and the RCDW approximation (Glass et al 1994), showed that for relativistic 
energies capture with spin flip is more probable than capture without spin 
flip. From the analysis above, we see that this is not the case in the extreme 
relativistic regime. 
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3 Conclusions 



In the Bates approximation logarithmic effects arise from scattering at the 
on-shell binary peaks. Our estimate of the corresponding terms in second- 
order perturbation theory confirm that there are indeed logarithmic terms, 
and the magnitude of these terms is larger than those obtained from the 
nonorthogonal correction introduced by Moiseiwitsch and Stockman (1980). 
This supports the claims made by Decker (1990) and Jakubassa- Amundsen 
and Amundsen (1985) in highlighting the cause of this effect and, more im- 
portantly estimates the coefficient of the logarithmic term. The total cross 
section including transitions involving spin flip is given by the formula: 

^robk 2 (7) ~ 2 5 ira u Z*Z 5 P (Z T + Z P ) 2 7 - 1 (ln 7 ) 2 (25) 

Interestingly, the logarithmic correction is not principally the result of retar- 
dation as previously supposed (Jakubassa-Amundsen and Amundsen 1985) 
but rather due to propagator singularities. In other words, it is a principally 
a kinematic feature. Furthermore the correction will only be significant for 
extreme energies and thus is primarily of academic interest. Lastly we note 
that capture without change of spin is the dominant transition at extremely 
high energies. 
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